Resonant nonlinear quantum transport for a periodically kicked Bose condensate 
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Our realistic numerical results show that the fundamental and higher-order quantum resonances of 
the (5-kicked rotor are observable in state-of-the-art experiments with a Bose condensate in a shallow 
harmonic trap, kicked by a spatially periodic optical lattice. For stronger confinement, interaction- 
induced destruction of the resonant motion of the kicked harmonic oscillator is predicted. 
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Dynamical systems that exchange significant energy 
with an external driving (time-dependent) field are 
paradigmatic objects for the study of complex evolutions 
with only a few degrees of freedom. In contrast to au- 
tonomous systems, where chaotic behavior can originate 
from many-body interactions, the complexity in driven 
systems arises from and can be controlled by the exter- 
nal field. An experimental setup where both types of 
complexity - many-body dynamics and external drive - 
are present is realizable and to a large degree controllable 
with state-of-the-art atom optical systems P, . Good 
control over nontrivial dynamics is the necessarytool for 
manipulating quantum states in a desired way Con- 
trolled coherent evolution is the necessary ingredient for 
quantum computing schemes, and, in particular, for al- 
gorithms which are based on a fast spreading over the 
Hilbert space of interest. Recently, ballistic expansion 
(i.e., linearly increasing momenta with time) was pro- 
osed to realize such quantum random walk algorithms 
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In this Letter, we answer the question of whether bal- 
listic resonant quantum transport can be realized with a 
periodically kicked Bose condensate. Cold dilute atomic 
gases have so far been used to realize many features of 
quantum chaos, such as dynamical localization j^j or dy- 
namical tunneling 0. To implement a fast spreading 
in momentum space, one can use quantum accelerator 
modes found recently 0, or the standard (5- kicked rotor 
(KR) which shows ballistic motion at the so-called quan- 
tum resonances "s^. The latter, however, are hard to 
realize if the initial conditions cannot be optimally con- 
trolled 9]. The preparation of a Bose condensate within 
harmonic traps offers very well-defined initial momenta, 
necessary for observing ballistic motion for a substantial 
number of kicks. 

We use the time-dependent Gross-Pitaevskii equation 
(GPE) [13 to describe a Bose-Einstein condensate in an 
harmonic confinement which is subject to a temporally 
and spatially periodic optical potential, created by a far 
detuned optical lattice. If the external potentials are not 
too strong, the GPE provides a good approximation of 
experiments with dilute Bose gases 1^- The GPE 



which we numerically integrate has the following form: 

+ ^ + ^ + Vo cos{2kLx) 

X Et=-oo F{t - mT) + gN\^{f, t)\^] ^(r, t) , (1) 

with ~ y"^ + z^- ilj{r,t) represents the condensate 
wave function, and M is the atomic mass. The nonlin- 
ear coupling constant is given hy g — ATrti^a/M, N is 
the number of atoms in the condensate, a the s-wave 
scattering length, and is the wave vector of the laser 
creating the optical potential. In principle, an arbitrary 
pulse shape F(t) may be realized, but here we restrict 
ourselves to a situation where the laser is switched on at 
time instants separated by T, with maximum amplitude 
Vq and periodic pulse shape function F{t) of unit am- 
plitude and duration ST <C T. Our system is then the 
nonlinear analogue of noninteracting cold atom experi- 
ments 0, 01, which implement the KR model, and we 
can directly compare our results with the well-studied 
KR (for ujx ^0), or with the Kicked Harmonic Oscilla- 
tor (KHO) 01 (for non- vanishing ujx)- The commonly 
used dimensionless kick strength and kick period of the 
KR are k = {Vq/H) dt F{t) and r = STEr/H, with 
the recoil energy = {hkLf /2M 0. 

The GPE was numerically integrated using a finite dif- 
ference propagator, adapted by a predictor-corrector loop 
to reliably evaluate the nonlinear interaction 14]. The 
external, time-dependent potential makes the integration 
challenging, in both time and computer memory. Typical 
integration times range from a few hours for a simplified 
ID version of Eq. (1), to several weeks for the full 3D 
problem. For the ID model, the motion is confined to the 
longitudinal {x) direction, and we use the renormalized 
nonlinear coupling parameter gio = 2hijj^a, assuming a 
radial trapping frequency ^ [Tsf. Experimentally, 
such a confinement is obtained using a cigar-shaped op- 
tical tube, as realized in the experiment of Moritz et al. 
[l^ . The initial state inserted into (1) is the relaxed con- 
densate wave function corresponding to the experimental 
ground state in a magnetic or optical trap. The ground 
state lies between the cases of a Gaussian, for 3 = 0, and 
the Thomas- Fermi limit, which is essentially an inverted 
parabola, for large nonlinearity Its characteristic 

width (Tp^ in momentum space is determined by the non- 
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linearity. Increasing the nonlinearity in Eq. (1) leads to 
a smaller width of the initial state (IQ] , 

The dynamics of the system described by (1) depends 
sensitively on the relative strength of the three poten- 
tials, i.e., on the control parameters Vq (for fixed pulse 
shape), uJx, and g. Since the system absorbs energy from 
the optical lattice, we must at all times compare the ki- 
netic energy with the (longitudinal) trap potential and 
the nonlinear term. If the latter two contributions are 
small, we can easily realize ballistic quantum motion, up 
to interaction times above which the trap potential is no 
longer negligible. On the other hand, we can tune the 
system to a situation where the trap crucially affects the 
dynamics, and we then have a realization of the KHO. 

In the following, we show that the fundamental as well 
as higher-order quantum resonances (QR) of the KR can 
be observed in an experiment using a Bose condensate, 
in the presence of a shallow harmonic confinement. In 
the linear KR, the QRs occur at specific kick periods 
T = Trs/r (s, r integer) ^. At the Talbot time T = Tt, 
the amplitudes of the wave function in momentum space 
exactly rephase between successive kicks for particular 
initial momenta (pi„it = 0) 0, Q . The result is a maxi- 
mal, i.e., perfectly phase-matched, absorption of energy 
from the kicks, leading to a ballistic spread of the wave 
packet [1,0. Only signatures of the QRs at T = Tt/2 
and T = Tt have been observed up to now in exper- 
iments with essentially non- interacting atoms, because 
the initial momenta of the atoms could not be sufhciently 
controlled 0]. 

Figure 1 presents our results for the fundamental QR 
of the KR at the Talbot time for Rb atoms, T = Tt = 
Trh/{2En) ^ 66.26 /is, with fci, ~ 8.1 x 10*^ m^^. Shown 
are the kinetic energy and the momentum distribution 
of the condensate along the longitudinal direction. The 
energy is computed from the momentum distribution, in- 
tegrated over the transverse directions, i.e.. 



E{K) = 
with P{px,K) = 



1 



dpxPlP{Px,K) 



dpydpz\ip{p, t = KT)\ 



(2) 



where K denotes the number of kicks. We present also 
the integrated spatial distribution along the transverse 
direction [in y, or equivalently z, because of the radial 
symmetry in Eq. (1)]: Py{K) = J dxdz\4>{f,t = KT)\^. 
As an example for higher-order QRs, which up to now 
have never been resolvable experimentally. Fig. 2 shows 
data for the resonance at T = Tt/4. For all the pa- 
rameters studied, the ballistic quantum transport, with 
a quadratic growth E{K) oc , is clearly visible. 

Deviations from the idealized ballistic motion arise be- 
cause of the contributions of the trap potential and the 
nonlinearity. At long times, the expanding condensate 
feels the harmonic trap potential, and further accelera- 
tion is hindered by the trap. This effect is negligible for 
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FIG. 1: (color online), (a) Kinetic energy (in units of 8-Er) 
YS K, &t T = Tt, k = O.l (pulse width 5T ~ 500 ns, rise 
time 70 ns; Vq/^Er ~ 8), = 10 Hz, ujJ2-k = 100 Hz, 

and A'^ — 10* (circles), = 5 x 10* (crosses). The linear 
KR evolution with pi„it = is shown by the dashed line, 
(b) Momentum distribution (2) after K = lh, N ^ 10* (solid 
line) , N = 5x 10* (crosses) . (c) The longitudinal (dotted line: 
K = 0; shaded area: K = 15) and (d) the transverse (dotted 
line: K = 0; squares: 7; plusses: 15) spatial distributions for 
A*' = 10* (thinner distributions) and A*' = 5 x 10* (broader 
distributions). The transverse dynamics is frozen due to the 
small effective nonlinearity. 



vanishing lu^. Equating the longitudinal trap potential 
and the kinetic energy, and using p^ ^ TrkKhk^ at T = 
Tt d, 01 , we estimate the kick number above which the 
trap dominates the dynamics as K ^ 2kLy/h/Mu>x — 55 
for the parameters of Fig.l. More crucial is the small but 
finite initial condensate momentum spread (cTp^ <C 2hk]^ 
P,0]), which, after a characteristic time t* oc l /gri , leads 
to a linear increase of the energy E[K) oc K Il7l|. This 
crossover sets in at if ~ 10 in Fig. 1 (a). For small Up^, 
the condensate ground state extends over many lattice 
sites of the kick potential, which in turn makes a smaller 
trap potential necessary for ideal ballistic motion. The 
interplay between these situations is illustrated in Fig. 3, 
where we systematically vary the number of atoms in the 
condensate. We stress that the effect of the nonlinearity 
manifests itself indirectly via its influence on the initial 
state, while the nonlinear interaction is negligible dur- 
ing the kick evolution. This finding is quite surprising, 
remembering that the QRs correspond to exact phase re- 
vivals between kicks. On short time scales, however, the 
perturbation induced by the nonlinearity cannot accu- 
mulate a large enough dephasing, because even for the 
small kick strength /c < 0.5 it is at least 1 order of mag- 
nitude smaller than the kinetic energy. Our results are 
consistent with a simplified ID model analysis of the QR , 
in the presence of a small nonlinear perturbation |l8l |. 
This analysis, however, could not account for the exact 
nonlinear wave packet evolution including the harmonic 
confinement. 
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FIG. 2: (color online). Kinetic energies ai T = Tt /4:, same 
trap and pulse shape as in Fig. 1: (a) ID, A'^ = 10^, k — 
0.1 (red diamonds) and k — 0.09 (blue stars), and (b) k — 
0.45 (Vo/En =i 42), = 10* (ID, solid line; 3D red stars), 
and 3D fc = 0.48 {Vo/En ~ 45), iV = 10^ (squares), (c) 
Momentum distribution (2) for 3D with k = 0.45 after K — 
15. (d) A scan of the kinetic energy vs T for the ID case from 
(b), after K = 10 (circles) and 7^ = 20 (red crosses), with 
the corresponding data for the linear KR (dotted and dashed 
lines) . The resonance at T = Tt /4 manifests itself clearly, and 
a tiny peak of another higher-order resonance is marked by the 
arrow. The inset in (a) confirms the correspondence between 
our method (circles) and a ID fast Fourier transform evolution 
of the linear KR (solid line) [a wave packet (circles) or an 
incoherent ensemble of plane waves (solid line) was evolved 
respectively, with initial Gaussian momentum distribution - 
(Tp^ = 0.026 hkL - k = 0.09]. For comparison, the linear KR 
(Pinit = 0) in (a) for A; = 0.1 (solid line), k = 0.09 (dotted 
line), and (b) k = 0.48, k = 0.45 (dotted fine). 



We finally note that for a large number of atoms 
A'' > 5 X 10^ the results of our ID model and the full 3D 
computations differ [c.f. Fig. 3(a)] because initially the 
condensate substantially expands along the transverse di- 
rections. This is crucial for the realization of ballistic 
quantum transport, since for the same number of atoms 
in the condensate, the trap has less effect in the 3D as 
compared with the ID case, but, on the other hand, tTp^ 
in the 3D case is slightly larger. For the parameters of 
Fig. 3(a), the larger cTp^ has a negligible influence, while 
the smaller initial spatial extension allows the ballistic 
motion to survive longer in the 3D than in the ID case. 

The nonlinearity manifests itself if we scan the kick 
period over the fundamental QR and plot the kinetic 
energy at fixed K. The resonance peak shows a slight 
asymmetry, which does not occur in the linear KR. The 
asymmetry decreases when (i) reducing the longitudinal 
confinement [solid line as compared with circles in Fig. 
4(a)], or (ii) evolving the initial condensate state with- 
out the nonlinear term [pyramids in Fig. 4(a)]. Again 
the trap potential hides the influence of the nonlinear- 
ity, but the observed asymmetry originates from both 
perturbations of the usual KR. Any such perturbation is 
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FIG. 3: (color online). Scan of the number of atoms (trap 
as in Fig.l if not stated otherwise): (a) for T — Tt, and 
ID iV = 5 X 10^ (crosses), N = 10* (circles) TV = 5 x 10* 
(diamonds), A'^ = 10^ (plusses); 3D in the inset for A'' = 
10^ (sofid fine), A^ = 5 x 10^ (crosses), A^ = 10* (circles) 
A' = 5 X 10* (diamonds). The dot-dashed line shows the 
analytical resuh E{K) = k'^K'^/A [1 inear KR, Pi„it = 0], and 
the dotted line ID data for LOxfiTT = 5 Hz, and N = 5 x 10*. 
(b) T = Tt/4:, fc = 0.45, ID for A^ = 10^ (stars), A^ = 10^ 
(solid line), A^ = 3 x 10^ (dashed), A^ = 10* (circles), A^ = 
5 X 10* (diamonds); 3D in the inset for A" = 10^ (solid line), 
Af = 5 X 10^ (crosses), A^ = 10* (circles). 
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FIG. 4; (color online), (a) Scan of the kinetic energy vs T 
around T — Tt, parameters as in Fig. 1, in comparison with 
the linear KR (blue dots; ap^ — 0.01 hkL), after K — 14; 
ID {N = 10*: circles; Lj^/27r = 5 Hz: solid line; A' = in 
kick evolution; pyramids), and 3D (A' = 10*: green squares; 
N — 5x 10*: red crosses), (b) Kinetic energy of 3D KHO with 
T = 250 Ais, uj^/2-n- = 1 kHz, u;,/27r = 5 kHz, and k = 0.1 for 
N ^0 (solid line) and TV = 100 (stars). 



expected to introduce an asymmetric peak shape, follow- 
ing a semiclassical analysis speciflcally developed for the 
description of decoherence by spontaneous emission |l7| . 

We can also tune our system described by the GPE (1) 
to the limit, in which the nonlinear interaction term be- 
comes very important and even dominant for small kick 
strength. Experimentally such a situation can be realized 
either close to a Feshbach resonance (where the scatter- 
ing length for two-body atomic collisions can increase 
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by orders of magnitude [lOj) or by simply increasing the 
strength of the harmonic confinement. The latter sit- 
uation has been analysed as a nonlinear generalization 
of the KHO in [llj, |l9j . The analysis showed clear signa- 
tures of the nonlinearity in the resonant transport regime 
of the KHO, i.e., for Tuj^ = 27r/r with r = 3,4,6, where 
the quantum transport is enhanced with respect to the 
classical one, much in the same way as at the QRs of the 
KR 01 . In the KR studied above, the total energy - 
corresponding to the chemical potential in the station- 
ary case |loj - is almost entirely given by the kinetic 
energy obtained from the kicks. In the KHO, the energy 
distribution is very different: here the kinetic, the poten- 
tial, and the self energy (given by the nonlinear term) 
are of the same order of magnitude, which leads to a 
nonlinearity-induced redistribution of energy also to the 
transverse degrees or freedom. 

Our numerical technique can be used to generalize the 
previous studies of the nonlinear KHO to the full 3D 
GPE, properly including the coupling of the transverse 
dimensions which is necessary for comparison with real- 
life experiments. Results on the KHO are presented in 
Fig. 4(b), where we observe the destruction of the reso- 
nant motion at Tlu^ = tt/2 already for a small number of 
atoms N — 100, and after just a few kicks. In contrast 
to what was done in [ill ITflj ] , the initial state need not 
be translated away from the classically stable origin in 
phase space. It suffices to use the relaxed ground state 
of the condensate to observe the impact of the now much 
stronger effective nonlinearity due to the strong harmonic 
confinement. 

We have tested our data in different ways, making sure 
that our numerical codes produce stable results. In the 
ID case and for giD = 0, we compared them either with 
analytical results for the linear KR , or with the much 
simpler evolution on a discrete grid, using a standard fast 
Fourier transform (FFT) [tR ll^ . The analytic growth 
rate of the kinetic energy for p^^^ = at T — Tt is shown 
in Fig. 1(a). For T = Tt/4, in the inset of Fig. 2(a), 
we compare both evolutions for up to if = 50 kicks. The 
FFT code does not take into account (i) the trap, (ii) 
the finite duration of the pulses, and (iii) the coherent 
evolution of the wave packet, which explains the tiny de- 
viations for K > 40. The damping of the oscillations is 
due to the nonzero , an effect which was experimen- 
tally observed also for the anti-resonance 
where the motion of the linear KR is perfectly periodic 
only for p^^^ ~ 0. 

In summary, we presented the first (3-1-1) dimensional 
treatment of a Bose condensate driven by an external 
temporally periodic force, which controls the dynamics 
of the system. Ballistic motion is shown to be realizable 



over a substantial number of kicks, even in the presence 
of a weak harmonic confinement. Within the framework 
of the 3D GPE, we have taken a first step towards the 
study of higher-dimensional chaos induced by the nonlin- 
ear coupling of the spatial dimensions, and future work 
will concentrate on situations where the transverse de- 
grees of freedom significantly contribute to the dynamics. 

This work was supported by MIUR, COFIN-2004, 
the Humboldt Foundation (Feodor-Lynen Program) , the 
Scuola di Dottorato G. Gahlei, and the ESF (CATS). 



[1] L. Deng et ai, Phys. Rev. Lett. 83, 5407 (1999). 

[2] G.J. Duffy, A.S. Mellish, K.J. Challis, and A.C. Wilson, 
Phys. Rev. A 70, 041602(R) (2004). 

[3] S. Potting, M. Cramer, and P. Meystre, Phys. Rev. A 
64, 063613 (2001); J. Gong and P. Brumer, Phys. Rev. 
Lett. 88, 203001 (2002). 

[4] See, e.g., D. K. Wojcik and J. R. Dorfman, Phys. Rev. 
Lett. 90, 230602 (2003); T. A. Brun, H. A. Carteret, and 
A. Ambainis, idid. 91, 130602 (2003). 

[5] M.G. Raizen, Adv. At. Mol. Opt. Phys. 41, 43 (1999). 

[6] D.A. Steck, W.H. Oskay, and M.G. Raizen, Science 293, 
274 (2001); W.K. Hensinger et al, Nature 412, 52 (2001). 

[7] M.K. Oberthaler et ai, Phys. Rev. Lett. 83, 4447 (1999). 

[8] P.M. Izrailev, Phys. Rep. 196, 299 (1990). 

[9] W.H. Oskay et ai, Opt. Comm. 179, 137 (2000); M.B. 
d'Arcy et ai, Phys. Rev. E 69, 027201 (2004); G. Duffy 
et ai, ibid. 70, 056206 (2004); S. Wimberger et ai, Phys. 
Rev. A (to be published), e-print physics/0502061 
[10] C.J. Pethick and H. Smith, Bose-Einstein Condensa- 
tion in Dilute Gases, (Cambridge University Press, Cam- 
bridge, 2002); L. Pitaevskii and S. Stringari, Bose- 
Einstein Condensation, (Oxford University Press, Ox- 
ford, 2003). 

[11] S.A. Gardiner et ai, Phys. Rev. A 62, 023612 (2000). 
[12] C. Zhang, J. Liu, M.G. Raizen, and Q. Niu, Phys. Rev. 

Lett. 92, 054101 (2004); B. Mieck and R. Graham, J. 

Phys. A 37, L581 (2004). 
[13] G.M. Zaslavsky et al., Weak chaos and quasi-regular 

patterns (Cambridge University Press, 1992); F. Bor- 

gonovi and L. Rebuzzini, Phys. Rev. E 52, 2302 (1995); 

A.R.R. Carvalho and A. Buchleitner, Phys. Rev. Lett. 

93, 204101 (2004). 
[14] E. Cerboneschi et al, Phys. Lett. A 249, 495 (1998). 
[15] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[16] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 

Rev. Lett. 91, 250402 (2003). 
[17] S. Wimberger, 1. Guarneri, and S. Fishman, Nonlinearity 

16, 1381 (2003); Phys. Rev. Lett. 92, 084102 (2004). 
[18] L. Rebuzzini, S. Wimberger, and R. Artuso, Phys. Rev. E 

71, 036220 (2005). 
[19] R. Artuso and L. Rebuzzini, Phys. Rev. E 66, 017203 

(2002). 



